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Abstract. Let H = A + J^ #a=2 V a be a 3-body Hamiltonian, H a the sub- 
system Hamiltonians, A > the Laplacian of the Euclidean metric g on 
Xq = IK™, V a real- valued. Buslaev and Merkurev ji b| have shown that, 
if the pair potentials decay sufficiently fast, for </> £ C£°(K), the operator 
<j>(H) — 4>(Hq) — '}2^ a —2(4>(H a ) — <P[Ho)) is trace class. Hence, one can define 
a modified spectral shift function cr, as a distribution on R, by taking its trace. 
In this paper we show that if V a are Schwartz, then cr is in fact C°° away 
from the thresholds, and obtain its high energy asymptotics. In addition, we 
generalize this result to TV-body scattering, TV arbitrary. 



1. Introduction 

Let H — A + X)#a=2 be a 3-body Hamiltonian, H a the subsystem Hamil- 
tonians, A > the Laplacian of the Euclidean metric g on X = R™, V a real- 
valued. Buslaev and Merkurev ||, ||] have shown that, if the pair potentials decay 
sufficiently fast, and under a spectral assumption, for € C£°(R), the operator 
4>(H) — (j>(Ho) — X)#a=2(^(-^a) — <j)(Ho)) is trace class. Indeed, this follows from the 



adoptation of the Helffer-Sjostrand functional calculus |lq| , see 1 30 , even without 
the spectral assumption. Hence, one can define a modified spectral shift function, 
which is really a distribution on R>, by 

(1.1) cx{4>) = tv{ct>{H) - 4>{H ) - Y, W H a) - 4>{B ))). 

#a=2 

In the two-body setting, a is often denoted by The purpose of this paper is to 
show that a is in fact C°° away from the thresholds, and to obtain its high energy 
asymptotics. Namely, we prove the following theorem. 



Theorem. (See Theorem 5.i.) Suppose that the pair potentials V a are Schwartz 
(on X a ). Then the spectral shift function a is C°° on R \ A, where A is the set 
of thresholds and L 2 eigenvalues of H. Moreover, a is a symbol outside a compact 
set, and it has a full asymptotic expansion as X — * +oo; 



a(A)~^At-3-i Cj . 
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In addition, 



C ° = C °EE / V a V b dg, 

a h=£n J X <> 



a b^a 

where Cq = jg(n — 2)(n — 4)(27r)~™ vol(§™ -1 ) depends only on n — dim Xq , and dg 
is the Riemannian density of the metric g. 

Buslaev and Merkurev had shown previously [pj that a is given by a continuous 
function, under the assumption that is not an eigenvalue or resonance of any two- 
body subsystem, and that the bottom of the essential spectrum is not an eigenvalue 
of the whole Hamiltonian. We remark that the leading term has an additional 
A -1 decay as compared to the corresponding formulae in two-body scattering, see 
HI HI m f§ . This is due to the fact that 4>{H) - 4>(H ) - E #a=2 0(#a) - 
(f>(H )) is lower order in a high-energy sense than <j>(H) — <p(H ) is in the two-body 
setting. Correspondingly, if we include a three-particle interaction V\ in H, i.e. 
allow j^a = 1, the leading term will have its usual order A^~ 2 . 

In fact, this is not the only modified spectral shift function one can consider. For 
example, let Xa, #a > 2, be a partition of unity by C°° functions on the radial (or 
geodesic) compactification Xo of the configuration space Xo that correspond to the 
collision plane structure. Then <j>(H) — J2# a >2 ^(-Ha)Xa ls & ^ so trace class, and one 
can define a modified spectral shift function by taking its trace. This spectral shift 
function will depend on the choice of the partition of unity, and it may be argued 
that it is less natural than the one adopted above. This trace is also smooth away 
from the thresholds and its high energy asymptotics can also be calculated. In fact, 
this expression generalizes without changes to many-body scattering with arbitrar- 
ily many particles, and our proof shows its smoothness. Even the Buslaev-Merkurev 
expression can be adopted to X-body scattering, N arbitrary, with combinatorial 
complexity being the only additonal issue. The corresponding result is stated in 



Theorem 5.3 in the last section. 

The main reason why these statements hold is that the corresponding expres- 
sions for the spectral measure, or for the high-energy cutoff outgoing (or incoming) 
resolvents tp(H)(H — (A + iO))^ 1 , V G C£°(R), already make sense pointwise in A. 
While these operators are not necessarily trace class, the trace makes sense as an 
oscillatory integral since it can be regarded as a pairing with the delta distribu- 
tion associated to the diagonal. Indeed, the latter only requires that the scattering 
wave front set of the kernel of the the operator whose trace we intend to take, 
is disjoint from the conormal bundle of the diagonal lifted to the b-double space 
[Xo x Xo;dXo x dXo]. In view of the propagation of singularities in many-body 
scattering, applied to the resolvent kernel, this is automatically satisfied provided 
that the kernel of the operator combination that we consider is essentially a general- 
ized eigenfunction of some many-body Hamiltonian microlocally near the conormal 
bundle of the diagonal. 

Our methods are thus an adaptation of traditional microlocal analysis to many- 
body scattering. In fact, we do not need the whole program initiated by Melrose 



21 1, see fl31| , since the trace is a rather simple object. Thus, a moderate strength- 
ening of propagation results in the dual of the radial variable essentially suffices. 
These have their origins in the Mourre estimate [ p4[ , and became partially microlo- 
cal estimates in the work of Gerard, Isozaki, Skibsted || [lCj and the work of the 
second author |B2| ; see also |13] for a discussion of these in the geometric setting. 
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The weak high energy asymptotics, i.e. that of 

tr(<Xff/A) - <f>(H Q /\) - i<t>{Ha/\) - HHo/X))), 

#a=2 

4> G (M) fixed, can be derived from the semiclassical functional calculus of Helffer- 
Robert |^7|, [l4| . Symbol estimates for a follow from combining microlocal versions 
of high energy estimates for R(\ + iO), |17| , with our method of proving smoothness 
of a (see also j33|). The perturbation series expansion for B,(\ + iO) then gives the 
full asymptotics. 

The assumption that the V a are Schwartz is not optimal. The distributional 
trace, a, is defined if the V a are symbols of order < — n (on X a ). Below we also 
obtain partial results for these potentials and explain the finer tools one needs to 
extend the Theorem to this setting. Specifically, this would entail obtaining, at 
least microlocally near the conormal bundle of the diagonal at infinity, oscillatory 
integral estimates in place of the wave front set estimates. 

We also mention that under stronger assumptions as in M , such as analyticity of 
the pair potentials in a cone near infinity, one can use complex scaling, or distortion 
analyticity, to continue (H — A) -1 across the spectrum. Since the trace is invariant 
under the scaling, and the scaled operators are elliptic near the real axis, we deduce 
that the trace is real analytic. 

While there is quite a bit of machinery in the microlocal analysis of this paper, 
due to the complexities of many-body scattering, given the microlocal estimates, the 
arguments proving the Theorem are essentially the same as in two-body scattering. 
Hence we invite the reader who is familiar with the scattering microlocal analysis 
introduced by Melrose J2lj to read Section || pretending that one works in the 
two-body setting. Moreover, in the appendix we sketch the proof of a propagation 
theorem that is weaker than the results of , but suffices for the purposes of this 
paper. The sketch is intended to stand on its own, and to explain the basic ideas 
in microlocalizing propagation estimates. 

The structure of the paper is as follows. In Section |^ we recall the standard 
many-body notation and make some preliminary remarks. In the following section 
we remind the reader of the microlocal structure of many-body scattering. In 
Section |J we prove the smoothness of a in the three-body setting, and in Section |^ 
we obtain the high energy expansion of a. In the final section we state and prove 
the corresponding results in many-body scattering, and in the appendix we outline 
the proof of the simplified positive commutator estimate described above. 

A. V. is grateful to Maciej Zworski for numerous stimulating discussions on traces 
in the geometric two-body setting (as arising in the work of Christiansen jij). The 
authors also thank Leonid Friedlander, Richard Melrose, Didier Robert and Steve 
Zclditch for helpful comments. 

2. Notation and preliminaries 

Before we can state the precise definitions, we need to introduce some basic (and 
mostly standard) notation. We refer to |6| for a very detailed discussion of the setup 
and the basic results. We consider the Euclidean space K™, and let g be the standard 
Euclidean metric on it. We assume also that we are given a (finite) family X of 
linear subspaces X a , a £ I, of R" which is closed under intersections and includes 
the subspace X\ = {0} consisting of the origin, and the whole space Xq = M™. Let 
X a be the orthocomplement of X a . We write g a and g a for the induced metrics on 
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X a and X a respectively. We let ir a be the orthogonal projection to X a , ir a to X a . 
A many-body Hamiltonian is an operator of the form 

(2.1) H = A + j2^ a yv a ; 

here A is the positive Laplacian, Vq — 0, and the V a are real-valued functions in an 
appropriate class. 

There is a natural partial ordering on / induced by the ordering of X a by inclu- 
sion, i.e. b < a means that X b C X a , hence Xb D X a . A three-body Hamiltonian is 
a many-body Hamiltonian with / ^ {0, 1} such that for all a, b ^ {0, 1} with a ^ b, 
I a nl(, = {0} holds. 

Corresponding to each cluster a we introduce the cluster Hamiltonian H a as an 
operator on L 2 (X a ) given by 

(2.2) H^Axa+^Vt, 

b<a 

Ax<* being the Laplacian of the induced metric on X a . With 

(2.3) X a — {X b n X a : b< a}, 

H a is a many-body Hamiltonian with collision planes X a . We write 

(2.4) H a = H a ® ld Xa + Id X a ®A Xa =A Xo +Y J V a . 

b<a 

The L 2 eigenfunctions of H a , also called bound states, can be used to define the 
set of thresholds of H b . Namely, we let 

(2.5) A a = U b<a spec pp (H b ) 
be the set of thresholds of H a , and we also let 

(2.6) A' a = A a U spcc pp (H a ) = U b < a spec pp (H b ). 

Thus, £ A„ for a ^ and A a C (-co, 0]. It follows from the Mourre theory (see 
e.g. |2^]) that A Q is compact, countable, and spec pp (H a ) can only accumulate at 
A Q . We also let 

(2.7) A = A 1 Us P ec pp ( J ff). 

For future reference we also define the intercluster interactions as 

(2.8) I a = H-H a = J2 V b- 

One of the properties of many-body operators that we exploit below is that if 
Hl and Hr are many-body Hamiltonians on (Xo)l> (Xq)r respectively, then 

(2.9) P = P a = aH L <g) Id+(1 - a)l&®H R , < a < 1, 

is a many-body Hamiltonian on Mq = (Xq)l x (Xq)r. (We name the new space 
Mq and the operator P to avoid confusing notation.) The collision planes are given 
by {X a )L x (Xb)R, a,b £ I, though they have a rather special structure generated 
by {X a )L x (X )r and (X )l x (X a )R, a G /, by taking intersections of these. The 
thresholds A = A Q of P are of the form chAl + (1 — cx)Xr, where Xl and are 
thresholds of Hl and Hr respectively. 

Since we want to study the Schwartz kernel of the resolvent on a space, for us 
the relevant case is X a = (Xo)l = (A"o)_r (this is the space), and Hl = Hr, so the 
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kernel of (Hl — A)~ x is a distribution on M$ = [Xq)l x {Xq)r. The thresholds of 
P a vary with a, except that A S = implies that A S A Q for all a. On the 
other hand, given A ^ A^, there is a countable subset C of (0, 1) such that a C 
implies A ^ A Q . Indeed, for any pair (A^, Xr) 6 Af,x Ar with Al ^ Xr, there is a 
unique a € K such that A = aA^ + (1 — a)A^, and then we can take C to be the 
intersection of the set of such a's with (0, 1). 

In the next section we explain the full microlocal propagation picture proved in 



31 1 . However, we would like to emphasize here that apart from mild technical issues 



we only need an estimate of Gerard, Isozaki, Skibsted ]9|, |10| and the second author 
32 1, and its slightly strengthened version. Namely let B = ■ D w + D w ■ j^j, 
and suppose that F 1 ,F 2 € C£°(M), suppF x C (c, d), suppF 2 C (d, c"), c < c' < c", 
\c\, \c"\ sufficiently small, and ip <E C^°(R). The aforementioned papers show that 
Fi(B)tp(H)(H - (A + iO))- 1 F 2 (B) maps H°' s to H°' s ' for all s,s' when A i A. 
In the terminology of the next section, this holds because the principal symbol 
of B at infinity, which is denoted by — r there, is increasing along generalized 
broken bicharacteristics, and singularities propagate in the forward direction under 
(H- (X + iQ))- 1 . 

We only need the following strengthening of this statement (as well as its high 
energy version). Suppose that xijX2 € S°(Xo) axe conic cutoff functions (outside 
a compact subset of Xq>) with suppxi n suppx2 = 0- Then there exists e > 
(depending on the supports and on A ^ A) such that for Fi,F 2 £ C£° (R) , supp F\ C 
(c, d + e), suppF 2 C (c', c"), c < c' < c", \c\, \c"\ sufficiently small, and ip € C£°(R), 
XiFi(B)ip(H ){H- (\ + iO)y 1 F 2 (B) X 2 maps H°' s to iJ°' s ' for all s,s' when A ^ A. 
In the microlocal terminology, this means that the principal symbol of B at infinity 
has to increase by at least e > along generalized broken bicharacteristics leaving 
suppx2 (near infinity) and ending in suppxi. We remark that the smoothness 
order r of the Sobolev space H r ' s is left as only due to the non-ellipticity (in 
the usual sense) of B. Slightly different pseudo-differential cutoffs, such as those 
described in the next section, yield maps from H r ' s to H r ' s for all r, r' , s, s' . 

We also need the high energy version of both of these estimates, which can be 
proved similarly. Namely, first, that for all r, Sjr'js', 

(2.10) \\Fi(B)<P(H)(H - (A + iO))- 1 ^)!!^,^,^ < CX' 1 / 2 , A > 1, 

if F U F 2 e C C °°(R), suppf\ c (c,c'), su PP F 2 c (c',c"), c < d < c", \c\,\c"\ 
sufficiently small, and rp € C£°(R). Next, there exists e > (depending on the 
supports of xi and X2 as above), but not on A > 1) such that for Fi,F 2 £ C£°(R), 
suppi^i C (c, d + e), suppF 2 C (c',c"), c < d < c", |c|, \c"\ sufficiently small, and 
ip S C~(R), and 

(2.11) WxiF^BMH^H - (A + *0))" 1 ^ 2 (i?)x2|| e(ff o. s , H ov ) < CA^ 2 , A > 1. 



3. Microlocal analysis in many-body scattering 

We recall the microlocal tools introduced in the many-body setting in . Thus, 
we compactify Xq as in pjj by letting 

(3.1) X = X = SI 

to be the radial compactification of Xq (also called the geodesic compactification) 
to a closed hemisphere, i.e. a ball. Identifying X Q with R n (and using g) we write 
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the compactification map RC : R™ — > S™ given by 

(3.2) RCH = (l/(l + |w| 2 ) 1/2 ,w/(l + |w| 2 ) 1/2 )eS^cR n+1 , 

One can view the use of the compactification as using inverted polar coordinates 
(r _1 ,w) € [0, f) x §™ _1 near infinity, i.e. working near 'r _1 = 0', where w = ru> in 
usual polar coordinates. We write the coordinates on R n = X a © X a as (w a , w a ). 
We let 

(3.3) X a = cl(RC(X a )), c a = x a n dS", 

Hence, C a is a sphere of dimension n a — 1 where rt Q = dimX a . We also let 

(3.4) C = {C a : a£ I}. 

Thus, C = dSl = S n -\ and a < b if and only if C h C C a . We also define the 
regular part of C a as 

(3.5) C ai r Cg — C a \ Ub^ a Ch — C a \ [_J{Cfc : Cb C C a } 

Since throughout this paper we work in the Euclidean setting, where the notation 
X, X a , etc., has been used for the (non-compact) vector spaces, we always use a 
bar, as in X, X a , etc., to denote the corresponding compact spaces. We write H r 
for the standard Sobolev space, and H r ' s for the weighted Sobolev space (w)~ s H r . 

In |]30| a pseudo-differential operator calculus was constructed on many-body 
spaces as above (M. n ,X); indeed, it was defined for appropriate geometric general- 
izations of these spaces. Operators of multi-order (k,l) are denoted by ^> s ' c (Xq;C); 
the notation only refers to C rather than X since this is the only pertinent infor- 
mation in the geometric setting. We briefly recall the definition of the many-body 
pseudo-differential calculus via the quantization of symbols. Below [Xo;C] is the 
blow-up of Xq at C; in particular, its interior is diffeomorphic to that of Xq, i.e. 
to Xq = R n , and there is a smooth blow-down map [-Xo;C] — > Xq, so every C°° 
function on X is C°° on [X ;C]. Thus, the polyhomogeneous space (Xo,C) is 
the following. We identify int(Xo) and int([X ;C]) with R™ as usual (via RC -1 ), 
suppose that a e C°°(R™ x R£) is in fact of the form 

(3.6) agp^r([l„;C]xl*), 

where pg and poa are defining functions of the first and second factors, X and Xq, 
respectively, so they can be taken as (w)" 1 and (O" 1 respectively. Let A = ^(a) 
denote the left quantization of a: 

(3.7) Au{w) = (27r)-™ J e l(w ~ w ' H a(w, Qu(w') dw' d£, 

understood as an oscillatory integral. Then A £ 4'^,' (Xq,C). We could have 
equally well used other (right, Weyl, etc.) quantizations as well, and we could have 
also allowed a to depend on w' as well (with w' regarded as a smooth variable on 
[X Q ;C]). 

An advantage of using symbols 

(3.8) a S C°°{[X -C] W x [X ;C\ W > x {X*) € ) 

is that certain operators with amplitudes a — q(w,£)p(w' which can be consid- 
ered the operator product A = QP of the left quantization Q of q and the right 
quantization P of p, lie in the class ty™/ (Xq,C) even though the left quantization 
of q does not. As an example, for ipo £ C£°(R), we can write ipo(H) as the right 
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quantization of a symbol p as above that is in fact Schwartz in £. Then q does not 
have to lie C°°([X ;C] x JJ) to make fl3.8| ) hold, since due to the rapid decay of p 
in £, the £ — > oo behavior of q is (mostly) irrelevant. For example, we can allow 



q = f(w- f/M)> / G then a = qp satisfies (gj). Thus, A = A(/) defined 

by 

(3.9) 

Au(to) = (2tt)-" / e^^'Ha^, u/, dw' d£, o(to, u/, = q(w, €)p(v/, £), 

yields an element of ^g c °°' (Xq,C). The operator has very similar properties 
to f(B)ipo(H), B as in the previous section, and can be considered as a pseudo- 
differential replacement for f(B)tpo(H), for f(B)tpo(H) is not a ps.d.o. (due to the 
lack of ellipticity of B). (Note though that it follows from the arguments of [|l3) 
that MH)f(B)MH) G *s c °°' (Xo,C) for Vi,</>2 G C~(E).) 

Elements oi^(X ;C) are in particular bounded operators from H r ' s to i? r . 
(The different signs are due to the index convention introduced in the geomet- 
ric scattering setting in 0].) It was also shown in p(J that for A ^ spec(-ff), 



R(X) = (H — X) 1 G ^g 2 ' (^o,C), and the Helffer-Sjostrand argument was adopted 
to show that for <f> £ C£°(R), <f>(H) G ^7°°' (Xo,C), so it is smoothing in the stan- 
dard sense, but does not give any decay at infinity (see Section || for a sketch of 
this argument in the semiclassical setting). The construction of R(X) is local at 
8Xq (in fact, it is microlocal in a sense discussed below). Thus, if H = A + V, 
H' = A + V are many-body operators such that in a neighborhood of p S dXo, 
H — H' = V - V is in x k C°°(X ), k > 1, i.e. x(^ - V"') G x fe C°°(Z ) for 
some x G C°°(Xo), x(p) 7^ °> then #ff(A) - i?#'( A ) is in *s c 2 ' fc (^o,C) near p, so 
(Rh{X) — Rh'(X))x G ^/g c 2 ' fe (X ,C). This in turn implies, via the Helffer-Sjostrand 
construction, that for <fi G C~(K), </>(#) - tf>(H') is in *g c °°' fe (Xo,C) near p. 

Applying this in our setting, i.e. for 3-body Hamiltonians H with potential 
^ #a=2 V a , note that if the potentials V& are classical symbols of order —k on 
X b , then H — H a is in x k C°°(X ) away from U#b=2, b^aCb, and H a — H Q is in 
x k C°°(Xo) away from C a . Hence 4>(H) — 4>(H a ) is in (Xq,C) away from 

U#b=2, b^aCb, and <fi(H a ) — 4>{Hq) is in the same class away from C a . Now consider 
the expression <f>(H) — <p{H ) — ^2^ a=2 (<j>(H a ) ~ 4>{Ho)) near Cfe- Rewriting it as 
4>{H) - cj>(H b ) - E #a =2, a^( H «) ~ shows that it is in *g c °°' fe (Xo,C) away 

from U^t Q= 2, a =jtbC a , hence in particular near Ct,. Since & is arbitrary, we deduce that 
4>{H) - <t){H ) - E# o =2(0C ff a) - 4>{H )) G ^s^°' k (X ,C). Since pseudo-differential 
operators of order (to, I), m < — n, I > n, are easily seen to be of trace class (the 
kernel is in particular continuous and integrable along the diagonal), we deduce 
that (j>(H) — 4>(H ) — J2^n=2( ( l ) (^a) ~ ^{Ho)) is trace class indeed, as discussed by 
Buslaev and Merkurev H. This allows us to define a by taking its trace, i.e. by 

& 

We need more refined tools to analyze the smoothness of a. Namely, we need 
to analyze the Schwartz kernels of operators that do not lie in ^S.' (Xq,C), e.g. 
of (H — (A + i0)) _1 . This requires the use of some microlocal analysis, which we 
briefly recall below. 

The phase space in scattering theory is the cotangent bundle T*Xq which can be 
identified with X x (X )* . Again, following Melrose M , it is convenient to consider 
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its appropriate partial compactification, i.e. to consider it as a trivial vector bundle 
over Xq by compactifying the base. Hence, one defines the scattering cotangent 
bundle of X by sc T*Xq — X x Xq. We remark that the construction of SC T*X is 
completely natural and geometric, just like the following ones, see pl| . Note that 
(Xq)* can be identified with Xq via the metric g as we often do below. 

Recall from that in many-body scattering sc T*Xq is not the natural place 
for microlocal analysis for the very same reason that introduces the compressed 
cotangent bundle in the study of the wave equation on bounded domains; see the 
works of Melrose, Sjostrand |2^] and Lebeau |l9j on the wave equation in domains 
with smooth boundaries or corners, respectively. We can see what causes trouble 
from both the dynamical and the quantum point of view. Regarding dynamics, the 
issue is that only the external part of the momentum is preserved in a collision, the 
internal part is not; while from the quantum point of view the problem is that there 
is only partial commutativity in the algebra of the associated pseudo-differential 
operators, even to top order. Hence, one cannot expect to localize in arbitrary 
open subsets of sc Tq^Xq, i.e. to microlocalize fully. To rectify this, we replace 
the full bundle sc T£ arog A = C a , rcg X (X )* ~ C a ^ cg x X over C a , lcg C dX by 
sc T£ a T ^X a = C Q , rcg x X a , i.e. we consider 

(3.10) SC T*X = U a sc T£ aics X a . 

Over C ajlcg , there is a natural projection ir a : SC J£ X — > SC J£ X a , i.e. ir a : 
Ca.rcg x Aq — > Ca.rcg x X*, corresponding to the pull-back of one-forms; in the 
trivialization given by the metric it is induced by the orthogonal projection to X a 
in the fibers. By putting the 7r a together, we obtain a projection ir : SC T*^^X — > 

sc T*Xq. We put the topology induced by ir on sc T*Xq. As mentioned above, this 
definition is analogous to that of the compressed cotangent bundle in the study of 
the wave equation in domains with smooth boundaries or corners. 

We also recall from that the characteristic variety So(^) of A— A is simply the 
subset of sc Tqx Xq where g — A vanishes; g being the metric function. If Ai = {0}, 
the compressed characteristic set of H — A is 7r(£o(A)) C sc T*Ao. In general, all 
the bound states contribute to the characteristic variety. Thus, we let 

(3.11) E 6 (A) = {C = (2/6,6) G sc T* c X b : A - |6| 2 G spec pp H b } C sc T^X b] 

note that |6| 2 is the kinetic energy of a particle in a bound state of H b . If C a C C&, 
there is also a natural projection 7Tf> a : SC J£ X b — > sc y£ x a (in the metric 
trivialization we can use the orthogonal projection X b — > X a as above), and then 
we define the characteristic set of H — A to be 

(3.12) S(A) = U a S a (A), S a (A) = U CbD c^ba^b(X)) n sc T£ a ias X a , 
so S(A) C sc f*X a . 

The radial subset R+(X) U -R-(A) of S(A) plays a special role in many-body 
scattering, for points in this form stationary generalized broken bicharacteristics. 
Namely, that at points (w/\w\,^) £ -R±(A), the Hamilton vector field of Ax b is 
radial (so there is no propagation tangentially to T*X b as discussed below), and 
simultaneously the particles may be in a bound state of H b , hence there is no 
propagation in normal directions either. More explicitly, let — r be the dual of the 
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radial Euclidean variable, so 

(3.13) T = -J^- 

Then 

R± (A) = U = (w, £«) e sc T5 a _ X a : 

' 4 ! 36, C a C C b , A - r(0 2 G spec^ff"), ±r(0 = |£ a |} 

are the incoming (+) and outgoing (— ) radial sets respectively. Note that if A is 
not a threshold of H, then r^Oon R+M U #_(A). The set R+(X) U i?-(A) is 
the propagation set of Sigal and Soffer [p9| ; this is the set where there is no real 
principal type propagation. 

Generalized broken bicharacteristics of a many-body Hamiltonian H — A were 
defined in Definition 2.1]. Here we do not recall the full definition here, but 
remind the reader that these are continuous maps 7 : I — > 11(A) , where J C K 
is an interval, that enjoy certain estimates with regard to Hamilton vector fields 
in the various subsystems. The main property for us is that the function r = 
— on sc T*Xo, is strictly decreasing along generalized broken bicharacteristics 
except at the radial points. Note that — r is the dual of the radial variable, and as 
mentioned in the introduction, it played a major role in previous works on many- 
body scattering. 

As mentioned in the introduction, 'singularities' (i.e. lack of decay at infinity) of 
u G S' are described by the many-body scattering wave front set, WFs c (u), which 
was introduced in pPf , and which describes u modulo Schwartz functions, similarly 
to how the usual wave front set describes distributions modulo smooth functions. 
Just as for the image of the bicharacteristics, sc T*Xq provides the natural setting 
in which WFs c is defined: WFs c (u) is a closed subset of SC T*X . The definition of 
WFs c (u) relies on the algebra of many-body scattering pseudo-differential opera- 
tors, $gj. (Xq,C). There are several possible definitions of WFs c , all of which agree 
for generalized eigenfunctions of H , but the one given in p0| ] that is modelled on 
the fibred-cusp wave front set of Mazzeo and Melrose [ p0[ enjoys many properties 
of the usual wave front set. 

We remark that in the two-body setting, when sc f*X = sc T* Xn X , WFs c is 
just the scattering wave front set WF SC introduced by Melrose, Q, which in turn 
is closely related to the usual wave front set via the Fourier transform. Thus, for 
€ sc T* Xq X , considered as dX Q x M™ = dX Q x R™, £ WF sc (u) means 

that there exists <f> € C°°(Xo) such that </>(w) ^ and J-((pu) is C°° near £. If we 
employed the usual conic terminology instead of the compactified one, we would 
think of as a conic cut-off function in the direction u>. Thus, WF SC at infinity is 
analogous to WF with the role of position and momentum reversed. The definition 
of WFsc(m) is more complicated, but if u = ip(H)v for some tp € C£°(R) (any other 
operator in ^g,, 00 ' (X,C) would do instead oftp(H)), then the following is a sufficient 
condition for (u, £ a ) e sc T£, a t X a , considered as C a;reg x X a , not to be in WFs c (u). 
Suppose that there exists (j> € C°°{X ), <t>(ui) ^ 0, and p € C%°(X a ), p(£ a ) ^ 0, and 
{{■K a )* p)F{(j)u) € S(W l ) = S{X ). Then (w,^) £ WF Sc (u). This wave front set 
also has complete analogues for the relative wave front set WF s j(u), i.e. when 
one is working modulo weighted Sobolev spaces H k ' 1 . In addition, although the 
definition of sc T*Xq depends on the metric g, there is a natural identification of 
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the compressed spaces corresponding to different metrics g, and the wave front set 
is defined independently of the choice of g then. 

As an example, consider Mq = Xq x Xq as the space, and let u G S' (Mq) be 
the Schwartz kernel of tp(H), H a many-body Hamiltonian on Xq, ip G C£°(R). As 
discussed above Mq is naturally a many-body space with collision planes I a x Ij, 
a, b 6 /. However, to accommodate the diagonal singularity, we introduce additional 
collision planes by adding the diagonal diag = {(w, w) : w G Xq} to this set, as 
well as the diagonals of X a : diag a = diagn(A a x X a ), a G /; the resulting set 
is denoted by M.. Recall that diag is the closure of diag in Mq, hence it is the 
b-diagonal as discussed in [3(J. Now the kernel K of ip(H) is conormal to 
<9diag a in the strong sense that it becomes C°° upon the blow-up of these (in the 
appropriate order); indeed, this is essentially the definition of ^s^°' (Xq;C). Note 
that if ipo G CJ°(R) is identically 1 on [inf supp ip, sup supp ip] , then ipo(P a )K = K, 
as follows directly from the functional calculus since H ® Id, ld®H commute, so 
the Fourier transform description of the wave front set given above is applicable. 

Since K is the kernel of a pseudo-differential operator, it is immediate that 
WFsc(-?Q C ^a sc T*-^g diag a , i.e. it lies over the boundary of the diagonal; 

indeed, the kernel is rapidly decreasing elsewhere by the definition of ^a^°' (Xo;C). 
We claim that for all a, WFs c (A") n SC T*— — diag„ lies in the zero section of 

SC T*^— diag a . To see this, let <j> G C°°(M ) be such that supp0 is disjoint from 

all collision planes other than those containing diag a . This means that <j> localizes 
near the regular, and away from the singular, part of <9diag a . The orthocomplement 
of diag a , with respect to <?l + gii, is A.i&g X a^ X a ®{{w,—w) : w G Xq}, and the 
corresponding orthogonal projections give coordinates (w a + w' a , w a + (w') a , w — w') 
(up to a factor of y2) • The corresponding dual variables are (up to another factor) 
(£a C + (C) a -, £ — Since <j)K is Schwartz vaw — w', and is symbolic inw + io' 
after some blow ups, in particular symbolic in w a + w' a , the Fourier transform of 
cj)K is Schwartz outside £ a + ^ = 0. Hence WF Sc (K) lies in £ +£' a = 0, i.e. in the 
zero section of SC T*— — diag„ as claimed. 

odiag aro(5 ou 

A similar argument proves that WFs c (<5diag) lies in the zero section of sc T*-r. — diag a . 
Indeed, one can see this directly from the definition, for any vector field tangent to 
diag annihilates <5di ag - 

Note that the orthogonal projection of the conormal bundle of diag to T* diag a 
is the zero section, and for this reason we call the zero section of SC T*— — diag„ 

adiag arcB oa 

the compressed conormal bundle of diag at diag a . We summarize these statements 
as a lemma. 

Lemma 3.1. The kernel K of ip{H), ip G C£°(R), satisfies WF Sc (K) is a subset 
of the zero section of SC T*— — diag a , i.e. of the compressed conormal bundle of 

the diagonal at diag a . The same holds for <5diag- In particular, r = on WFsc(K) 
and on WFs c (£ diag ). 

Remark 3.2. For example, consider A(F) G * Sc °° ,0 (M , dM) defined by (Q, F G 
C£°(R), and note that on its operator wave front set (a notion that we briefly recall 
in the Appendix; see |^, Section 5] for a detailed discussion) — r G suppF. The 
lemma immediately implies that for F G C^°(M) with ^ suppF, A{F)K G S. 
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The main property of wave front sets that we use is the following. If it, v g S' , 
WF Sc (u) is a compact subset of sc f*X , and WF Sc (u) n WF Sc (v) = 0, then the L 2 
pairing extends by continuity to define (u,v). For the relative wave front set one 
has stronger results. In particular, it suffices that WFg^(u) n WFj^ - («) = 0; the 
smoothness order k is irrelevant since WF§ c (m) is assumed to be compact. This 
property follows immediately from |?(| Proposition 5.4] and the definition of WFg c . 

The theorem on the propagation of singularities is the following. 



Theorem. (|3lJ Theorem 2.2}) Let u £ S'(R") ; A e R. Then 



(3.15) WF Sc {u)\WF Sc {(H -A)u) 

is a union of maximally extended generalized broken bicharacteristics of H — X in 
E(A)\WF Sc ((fl--A)«). 

A stronger version of this is valid for u = (H — (A + iO))" 1 /, defined if WFg c (/) 
satisfies appropriate non-incoming conditions (see Section 2] for a full discus- 
sion). Namely WFs c (/) only propagates forward along generalized broken bichar- 
acteristics. The non-incoming conditions are satisfied, for example, if r < on 
WFsc(/), and then we conclude that r < on WFs c (w). These statements are 
also proved by positive commutator estimates, and consequences can be easily re- 
stated as operator estimates. For example, consider A(F) € (Xq, C) defined 
by ([5j]). The fact that — r is increasing along generalized broken bicharacteristics 
then converts into the estimate of [?2j . Indeed, mapping properties of the wave 
front set, J|(| Proposition 5.4], imply that for A in a compact subset of R \ A, 
r, s,r' , s' arbitrary, 

(3.16) \\A(F)R(X + iO)A(F)\\ B{Hr ^ Hr ,, sl) < C, 

where F <E C£°(M), and in addition suppF C (c, d), suppF C (c',c"), c < d < 
c". |c|, \d\ sufficiently small. The fact that in addition — r is strictly increasing 
on any non-stationary generalized broken bicharacteristic implies the strengthened 
estimate that for all XiX € C oo (A > ) with suppx H suppx n dX = 0, and for A in 
a compact subset of R \ A, r,s,r' ,s' arbitrary, there exists e > such that 

(3.17) \\xA(F)R(X + iQ)A(P)x\\ B(m ,. iB r>,.> ) < C, 

where F S C^°(R), and in addition suppF C (c, d + e), supp# C (d , c"), c < d < 
c". |c|,|c'| sufficiently small. These are very similar to the estimates mentioned 
at the end of Section 0. Of course, these estimates can be simplified, but here we 
explicitly wanted to emphasize the connection between WFs c and the estimates. 
We outline the proof of ( |3.17 ) in the appendix, so that the reader can avoid the 



more complicated arguments leading to the stronger results of |5l| 



4. Smoothness of the trace 

Let K = K x , K a = (K a ) x denote the kernels of R(X + iO) = (H - (A + i0))"\ 
and R a {X + iO) = [H a — (A + iO)) -1 respectively, in Im A > (of course +i0 can be 
dropped if ImA > 0). Hence they are tempered distributions on X x X . Then 
K solves (Hl — X)K = <5di ag , {Hr — X)K — <5dia g , where Hl, Hr act on the left 
and right factors of Xo x Xq respectively, and <5di ag is the kernel of the identity 
operator, hence is the natural Dirac delta distribution associated to the diagonal. 
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(The second equation really contains the transpose of Hr — A, but since Hr is 
self-adjoint and real, this is just Hr.) Thus, for example, 

(4.1) (aH L + (1 - a)H R - \)K X = 5 diag , < a < 1, 
for ImA > 0. Hence, for ImA > 0, 

(4.2) K x = (aH L + (1 - a)H R - A)" 1 ^. 

As remarked in Section ^ P a = ai?^ + (l— o)Hr is a many-body Hamiltonian, so all 
the usual results on many-body Hamiltonians, such as propagation of singularities, 
are applicable. In particular, if A ^ A, we can arrange that A (fi A a by an appropriate 
choice of a. We let r = — (w •£ + «/• £')/((w, w')) as usual, with (w, w') denoting 
coordinates on Xq x Xq, and (£,£') the dual coordinates. Thus, — t is the dual of 
the radial variable r = \(w, w')\ with respect to the metric gL+9R- We thus deduce 
that (aHi, + (1 — o)Hr — (A + i0))^ 1 (5di ag makes sense since r = on WFs c (^diag) 
by Lemma |3.l| . Thus, 

(4.3) K x = (aH L + (1 - a)H R - (A + iO))" 1 ^ = (P a - (A + zO))- 1 ^. 

For technical reasons it is convenient to eliminate the diagonal singularity (which 
cancels from the spectral measure anyway). This can be done by considering 

(4.4) ij){H)(H- (A + zO))- 1 = $(H) 

instead of R(\ + iO), where i/j £ CJ°(R) is identically 1 near A. This has similar 
properties to R(X + iO), except that it does not have a diagonal singularity in the 
interior. In particular, its kernel, which we denote by K — K\, satisfies (Hl — 
X)K = (Hr — \) K = K^h)i K-4>(H) denoting the kernel of ij)(H). Hence, again 
using Lemma |3.l| , 

(4.5) Kx = (aH L + (1 - a)H R - (A + iO))" 1 ^^) = (P a - (A + iQ))- l K m) . 

As recalled in the previous section, r is decreasing under the forward generalized 
broken bicharacteristic relation. Then, by the propagation of singularities, namely 
that they only propagate in the forward direction under (P a — (A + iO))^ 1 , we 
deduce that r < on WFg c (^)- Similar statements hold for K a . Hence r < on 
WFscOT), 

(4.6) K' = K-k Q -J2(K a -K ). 

a 

We can write this explicitly by saying that if F £ C£°(R), supported in (— oo,0), 
then A(F)K £ S, where now A(F) is an operator on Mq = Xq x Xq, with prin- 
cipal symbol at infinity the 'same' as that of F(B)ip(H) (though the latter is not 
quite a ps.d.o.), B the radial vector field on X n x Xq, with principal symbol — r. 
Moreover, r is strictly negative everywhere on this wave front set except possi- 
bly at the compressed conormal bundle of the diagonal, i.e. at the zero section of 

We remark that while microlocally on WFs c (^), K is not Schwartz, it is of 
course still in a weighted Sobolev space, namely in i^ r > s ~ 1 for any s < — n/2 and 
for any r. This holds because K^^ H ) is in H 7%s for any s < —n/2 and for any r, 
and (P a — (A + iO)) -1 loses one order of decay. (This is the standard propagation 
phenomenon as for the wave equation, and the factor hr 1 appears in Q5.27 ) for the 
same reason.) 
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We claim that in fact r < tq < everywhere, hence in particular the compressed 
conormal bundle of the diagonal is disjoint from it. The key step is the following 
lemma. 

Lemma 4.1. Let b be a 2- cluster. ThenWF Sc (K-K b ) andWF S c(K a -K ), a^b, 
are disjoint from the zero section of SC T*— — diag,, and of SC T*— — diag n . 

d&agb.reg 9dla g ,re 9 1 

Proof. To see this, we first compute {P a — X)u, u = K — K b . Thus, 
(4.7) (P a - X)u = K^ (H) - K. HHb) + [a(I b ) L + (1 - a)(I b ) R ]K b . 



Thus, near 9diag b , the result is rapidly decreasing since the intercluster interactions 
{h)L, {h)R are such, and the same holds for the difference K^ih) ~ ^WiTO' It * s 
only here that we use V a <E S(X a ); see the comments of the following paragraphs 
for potentials that do not decay rapidly. Hence, in this region, WFs c (u) is a union 
of maximally extended generalized broken bicharacteristics of P a — A. That is, if ( 
is in SC T*— diagT U SC T*— diagn", and C G WF SC (u) , then there exists T > 

O dla g6,rog ° dla gO,rcg 

and a generalized broken bicharacteristic 7 : [— T, T] — > £(A) such that 7(0) = C 
and 7(t) S WFs c (u) for all t 6 [—T,T]. (The constant T depends on £, namely on 
how long it takes for bicharacteristics to leave SC T*- — diag,, \J SC T*— — diagv,.) 

9 dla gb,rog 9 dla g0,rcg 

But if £ is in addition in the zero section of SC T*— — diag b or of SC T*— — diag , 

c ,dla gb,rcg ddmg 0rEg 

then t(7(0)) = 0, and r( 7 (i)) > for t < 0, so j(t) <£ WF Sc (u) (as r < 

on WF Sc (u)) contradicting that ~/{t) S WF Sc (u) for all t G [—T,T]. Here we 

use again that A is not a threshold, for this assumption ensures that there are no 

stationary bicharacteristics at r = 0. Thus, the zero section of SC T*- — diag>. and 

ddiagf, rog 00 

of SC T*— — diag n , i.e. the compressed conormal bundle of the diagonal, is disjoint 

ddiag _ rog ™' 

from WFsc (u) . A similar argument shows that for b ^ a, WFsc(K a — Ko) is disjoint 
from the compressed conormal bundle of the diagonal in the same region. □ 

Combining the statements of the lemma shows that WFs c (^') has the same 
property for every b. Since b is arbitrary, this means that WFs c (^') is (globally) 
disjoint from the compressed conormal bundle of the diagonal, and r is negative on 
the former. (The existence of a strictly negative constant To such that r < tq on 
WFsc(^') follows by compactness.) In explicit terms, this means that there exists 
e > such that if F e C£°(R), supported in (-00, e), then A(F)K' e S. 

The trace of K' , defined as the push- forward of the restriction of K' to the diag- 
onal, thus makes sense by wave-front set considerations. Note that the trace is thus 
just the L 2 -pairing of K' with <5diag- Hence the only property of WFs c we need is 
that distributions with disjoint WFg c can be L 2 -paired, which we already discussed 
in the previous section. Again, explicitly this argument amounts to writing 

(4.8) 

(K',S diae ) = ((Id-^o(Pa))tf',<W> + (A(F)K',S dias ) + (K , ) (ik(P a )-A(F)*)5 atag ), 
F as in the previous paragraph, with in addition F = 1 near 0, and noting that 
both terms on the right hand side are defined. (Rather than using continuity from 
S in the style of Hormander Jn|, one could use continuity from Im A > here.) 

Note that this argument also shows that the trace of Xb{^{H)R{X + i0) — 
ip(H b )R b (X + i0)) is already defined as an oscillatory integral, or rather by wave 
front considerations. 
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Applying the same arguments for i?(A — iO), we deduce the same pointwise for 
the corresponding expressions involving the spectral measure, i.e. for 

(4.9) sp(A) - s Po (A) - 5> Pa (A) - sp (A)), 

a 

where 

(4.10) s P (A) = (2ni)- 1 ((H - (A + iO))- 1 - (H - (A - iO))" 1 ), 

and similarly with the other spectral measures. Since 4>(H) = J R 0(A) sp(A) dX, we 
conclude that the distribution a is given by the continuous function 

(4.11) a(X) = tr(sp(A) - sp (A) - $> Pa (A) - sp (A))). 

a 

The smoothness of this function follows directly by taking derivatives of the re- 
solvents with respect to A. Indeed the derivative of K with respect to A is given 

by 

(4.12) -^-Kx = (P a - (A + zO))- 2 ^^) = (P a - (A + zO))- 1 ^. 
Then t < still holds on WF Sc (^K x ), WF Sc (K x ). In addition, 

(4.13) ( Poi -\)JLk x = k x 

dX 

shows that 

(4.14) (P a - X)^(K X - (K b ) x ) =K\- {K b )\ + Wb)L + (1 - a){I b ) R ]{K b ) x . 



Since near <9diag b , WFs c (-Ka — (Kb)\) is in r < r < 0, we deduce as above that 
WFsc(^(^a — (Kb)\)) is also in r < r < near <9diag b , and then further that 
WFs c (-jtK') is in t < to < everywhere. This shows that ^ is continuous, and 
an iterative argument yields that a is C°° . We have thus proved the following 
proposition, which forms the first part of the Theorem stated in the introduction. 

Proposition 4.2. Suppose that H is a three-body Hamiltonian, and the pair po- 
tentials V a are Schwartz (on X a ). Then the spectral shift function a, defined by 
(0), is C°° onR\ A. 

If the potentials V a are not Schwartz on X a , rather they are symbols in S~ p (X a ), 
the first part of (5.27), namely K^,rjj) ~ ^V'(#&)> ^ s m H r ' s+P for any s < —n/2 



(and any r) near 9diag & , while [a(Ib)L + (1 — ct){h)R\Kb is in H r ' s ~ 1+P in the same 
region. This yields that u = K — K b is in jj r - s - 2 +P near the compressed conormal 
bundle of diag & (rather than Schwartz) . The pairing with <5diag thus makes sense if 
(s— 2+p)+s > 0, i.e. p > — 2s + 2. Since s > —n/2 is arbitrary, this means p > n + 2 
suffices. This result is non-optimal, p > n should be sufficient, but the improvement 
should correspond to more delicate (oscillatory integral type) behavior that we do 
not consider here since proving the better estimates is more complicated, especially 
if the set of thresholds is large! In addition, each derivative causes an additional 
order of decay to be lost, so to conclude that a is C k by this argument we need 
p > n + 2 + k. 

We explain briefly in terms of the two-body setting why more precise, oscillatory 
integral type, control is needed on the resolvent kernel near the conormal bundle of 
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the diagonal to obtain optimal results. In two-body scattering, the high-energy cut- 
off resolvent kernel K is of the form e l ^( w ~ w ) ( w _ u/)~( n_1 )/ 2 a, a smooth on the 
blown up space [Mo; <9diag] microlocally near the conormal bundle of the diagonal, 
and is also a smooth function of A > 0. Differentiating the kernel, in particular the 
exponential, with respect to A, gives an extra factor of growth (w — w ). However, 
this factor is bounded near the diagonal, so the pairing with <5di ag makes sense 
without imposing extra vanishing conditions on V. This argument shows that p > n 
suffices in the two-body setting to conclude that a is C°°; indeed, it immediately 
generalizes to the geometric setting of pl| , where it was proved by the use of trace 
formulae in 0], since one also knows the similar oscillatory behavior there |fr2"f . 
(This FlO-type analysis originated with J2^| in the geometric setting.) Such a 
detailed analysis is much harder in the many-body setting. 

Returning to the setting of Schwartz pair interactions, we remark that a perhaps 
better approach, which we do not pursue here, for proving smoothness, is to include 
A as a variable, i.e. to work with, say, the Schrodinger equation directly. Then the 
result of the pull-back to the diagonal in the spatial variables followed by the push- 
forward can be seen to be smooth by the wave front set calculus. 

We also emphasize that apart from technical issues, the whole argument can be 
rewritten in terms estimates involving operators such as F(B), B as in the previous 
sections. 



5. High energy asymptotics 

Below we change high energy problems into semiclassical problems, so we start 
by recalling some notation. First, H^' s = H^ s (M. n ) is the semiclassical weighted 
Sobolev space, (w)~ s HT(M. n ), where the norm on H£ is defined by 

(5-1) \\u\\ H r = (2nh)- n ^(0 r F h u\\ L 2 {Mn) , 

Th being the /i-Fourier transform, 

(5.2) (*fcu)(0= / e- iw ^ h u{w)dw. 



For r = 0, the norm in £2T(R") is the standard L 2 (R") norm; the norm on _ff^(R") 
is ||(1 + /i 2 A)w||i2( Rn ), etc. Semiclassical many-body pseudo-differential operators 
A = A(h) G ^>^/ h (Xo,C) are defined by the following analogue of (3.7): 

(5.3) (2nh)- n [ e t{w - w ' H/h a{w^;h)di 



where a S j3 00 m pgC 00 ([Xo;C] x Xq x [0, 1)), rapidly decreasing in the second factor. 
One can similarly modify (|3.9|), to define A — A(F; h) by 



(5.4) A(F; h)u(w) = (2nh)- n / e l( - w - w '>^ h a{w, w' , £; h)u(w') dw' d£, 



a(w , w' , £; h) = q(w,£-,h)p(w',£\h), this yields an element of ^g^' (Xq,C). One 
can also interpret these spaces as high energy spaces by letting h — A -1 / 2 ; then 
A — > oo corresponds to h — * 0. 

As an example, we compute the norm of kernel i^O) oiA = A(h) G ^gj^' (Xo,C) 
in H%"((X ) W x (X ) w >) for s < -n/2. Namely, takes the form 
(5.5) h- n A(w,(w-w')/h;h), A G C°°([X ;C] x X Q x [0, 1)), 
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A rapidly decreasing in the second factor. Thus, for s < —n/2 

\\K A{h)) f H0 „=h- 2n f ((w,w')) 2s \A(w,(w-w')/h;h)\ 2 dwdw' 



1. ~'2rl I {( nl . ,,J i\2.s| | ( ,,. , ,., ./],. /, ,|2 

- ,. , " Xo xXq 

(5.6) 

= h~ n / (w) 2s \A(w,W;h)\ 2 dwdW < C'h~ n , 

JXoxXa 

hence || Ka(H) \\h"- b — ChT n l 2 . In fact, a similar calculation, after inserting (Id +h 2 A(^ w w ^) r 
in front of A, gives that \\K A(h)\\H r ' a — Chr n l 2 for all r (and for all s < —n/2). A 
similar calculation also applies to 5diag, and gives for s < —n/2 that 

(5.7) Pdiaglk- ^ Ch- n ' 2 . 

We start our study of high energy asymptotics of a by stating an averaged 
version. 

Proposition 5.1. Suppose that H is a three-body Hamiltonian, V a are symbols of 
order -k on X a , k > n. Then for 4> G C^{R), 

(5.8) 

OO 

tr{4>(H/X) - <t>(H /\) - Y,(<t>(H a /X) - <t>{H Q /X))) - £ X^-^d,, A ^ +oo. 

Proof. We change the notation slightly and compute 

(5.9) tr(<f>(H/\') - <t>{H Q /\') - YM&JX) - 4>(H /X'))) 

a 

as A' — > +oo. As usual, we convert the high energy problem into a semiclassical 
one by letting h 2 = (A') -1 , and let 

(5.10) H(h) = h 2 H, 

so cf>(H/X') = <f>(H(h)), i.e. we are interested in the intersection of the spectrum of 
H(h) with a compact interval (namely supp^). Note that the high energy asymp- 
totics A' — > +oo corresponds to the semiclassical problem h — > 0. Note also that 
h 2 H = h 2 A + h 2 V, so semiclassically the potential vanishes two orders higher than 
the semiclassical Laplacian ft, 2 A. Let 4> £ C^°(C) be an almost analytic extension 
of (j) (so |<9 A <HA)| < CWIImAI^ for all N), and let R(X;h) = {H(h) - A)" 1 . Via the 
Cauchy formula, 

(5.11) <f>(H(h)) = — — [ dx4>(X)R(X; h) dX A dX, 

2m J 

the study of the functional calculus in a semiclassical setting reduces to that of the 
behavior of the resolvents away from the real axis, and the uniformity in h, up to 
the real axis in A. The main fact that makes the calculations uniform in h is that 
by self-adjointness, the L 2 operator norm of R(X;h) is bounded by |ImA| _1 , for 
this implies (via the parametrix construction) that all seminorms of i?(A; h) in the 
semiclassical version of *sc 2 '°(^o,C) (cf. (U)) are bounded by C\ Im A|~ fc for Re A 
in a compact set, | Im A| bounded from above, and for some C and k depending on 
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the seminorm. Hence it follows immediately that for <j> G C° 
(5.12) 

tr(<t>(H(h)) - 4>(H Q (h)) - Y,(<t>(Ha(h)) - cb(H (h)))) 



= tr (^-7^- J dx4>(\)(R{\: h) - R (X; h) - ^(i? a (A; h) - R (X; h))) dX A dX^J 

OO 

3=0 

While the integral is trace class, this may not be apparent since the integrand, 
which is in ^>g 2 ' k (Xo,C), is not such due to the diagonal singularity of the resolvents. 
However, it is easy to rewrite the integral so that the integrand becomes trace class 
as well. Namely, let ip(X) = (A - A o ) m 0(A), m > (dimX - 2)/2, A ^ R, so 
<fi(H(h)) = i?(Ao; h) m ip(H(h)), and write out the Cauchy integral representation of 
ip(H(h)) using an almost analytic continuation ip of ip. Thus, 

<P(H(h)) - 4>(H (h)) - ~ HHoW)) 

a 

(5.13) = / B^(X)(R(X : h) m R(X: h) - R Q (X ; h) m R Q {\; h) 



2m 



^2(R a (Xo; h) m R a {\- h) - i?o(A ; h) m R (\; h))) dX A dX. 



The integrand is now clearly trace class, in fact it is a semiclassical many-body 
scattering pseudo-differential operator of order (—2m — 2,p), and the asymptotics 
as h — > follows from the construction of R(X; h), etc. 

Our first remark is that since H is a differential operator (rather than a more 
general pseudo-differential operator), for all odd j, cij = 0, as in p6[ . This follows 
from the fact that these terms are given by integrals of odd functions of the mo- 
mentum £, hence these integrals vanish. Indeed, due to the additional h 2 vanishing 
in V , which implies similar vanishing for the difference of the various resolvents, the 
leading term is at least two orders higher, i.e. ao — (and we have already men- 
tioned that ai =0). This is exactly the same order of vanishing as for two-body 
scattering. In fact, more is true. 

To see this, we perform local calculations near Cf,, i.e. we replace the full trace 



by the trace of Xb times ( 5.13 ) where Xb G C°°(Xq) is supported away from 
the C a for a ^ b. Making sure that the x.b form a partition of unity, ( |5.12 ) 



becomes the sum of the local traces, so it suffices to show that the 02 term 
of each local trace vanishes. Thus, Xb(R(Xo', h) m R(X; h) — i?b(Ao; h) m Rb(X; h)), 
Xb{R a (Xo', h) m R a (X; h) — Rq(Xo; h) m Ro(X; h)), a ^ b, are already trace class. Using 
R(X; h) — Rb(X; h) = R(X; h)h 2 IbRb(X; h), a similar formula for R a (X; h) — Ro(X; h), 
that R(X; h), Rb(X; h) can be replaced by Ro(X; h) up to an error of h 2 , and the 
various operators can be commuted up to an error with an additional power of h, 
it is easy to see that the a-i term vanishes, where we also take into account that 
h = J2 a ^b^ a - (I n addition, 03 = since the odd coefficients vanish.) Changing 
back to the /i-independent notation finishes the proof. 

An alternative, but somewhat formal calculation (in the sense that it does not 
use the regularization procedure) proceeds as follows. Note first that the resolvents 
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can be rewritten as 

R(X; h) = R Q (X; h) - h)h 2 V a R (X; h) 



(5.14) 



Hence. 



(5.15) 



+ X R ( X '> h)h 2 V b R a {\; h)h 2 V a R (X; h), 

a b=jta 

R a (X; h) = Ro(X; h) - R a (A; h)h 2 V a R Q (X; h). 
R(X; h) - R {X; h) - £(i? a (A; h) - R {X: h)) 

a 

= X X R ^ h ) h2 V b R a (\ h)h 2 V a R (X; h). 



Thus, the combination of ( |5.15 ) and a regularization argument as in (5.12), though 
it is rather cumbersome to keep track of the regularizing factors, shows that there 



is an additional h 2 vanishing here, i.e. a<i — 0. Note that (5.15) is trace class when 
n = dimXo < 6, hence for two-dimensional particles, so in this case we do not need 
the regularization procedure. □ 

We proceed to find the leading coefficient, do = »4, without using the regulariza- 
tion procedure, which would complicate the calculation. This is justified directly 



for two-dimensional particles as mentioned above (the integrand of (5.12) is trace 



class); in general, we need to be more careful. To avoid cumbersome notation, we 
emphasize the case of dimX = 6, i.e. the particles are three-dimensional. First, 
one can write R(X; h) in a (finite) perturbation series: 



(5.16) 
R{X;h) 



R {X; h) - R {X; h)h 2 VR {X ; h) + R(X; h)h 2 VR {X; h)h 2 VR (X ; h); 



of course the last term can be expanded even further. Similarly, R a {X; h) can be 
expanded in a finite series. The sufficiently high order terms of the series (e.g. 
from the second term on if n = 6, i.e. if the particles are three-dimensional), when 
R(X; h), etc., are substituted into the right hand side (5.15), give trace class terms 
which also have at least additional vanishing in h compared to the leading term, so 
they do not contribute to do = 04. So for n = 6, 



(5.17) (5.12) = tr 



8\(f>(X) J2 J2 R °( X > h ) y bM^ h)V a Ro(X; h) dX A dX); 



we a priori know that the integral is trace class since it is the difference of trace 
class operators. In addition, the factors of V a , Vb, can be commuted to the front, 
and each commutator gives an extra h vanishing, and lowers the differential or- 
der of the pseudo-differential operator. In fact, the commutators take the form 
[R {X; h) , h 2 V a ] = -R (X;h)[h 2 A,h 2 V a ]Ro(X;h), so for n = 6 the corresponding 
terms are trace class with extra vanishing in h, so they do not contribute to do 
either. Thus, 



(5.18) 



= tr (~^ / ^( A ) X Y VbV a Ro(X; h) 3 dX A dX) 

1 Y Z7TI J — ' A — ' 
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As before, the integrand of Q5.18 ) is not trace class, though the integral is. Now 
using the Cauchy-Stokes formula, the integral can be written as 

(5.19) (Hl) = tr(|^ f 0(A)(^^HK)(i?o(A + iO;/ l ) 3 -i?o(A-iO;/ l ) 3 )dA), 

Now the integrand is not trace class, but its trace is well-defined as a pairing of 
its kernel with the delta distribution associated to the diagonal. We write out this 
pairing explicitly as the integral of the restriction of the kernel to the diagonal by 
writing Ro(X + iO; h) 3 — Ra(X — iO; h) 3 as multiplication by a differentiated delta 
distribution associated to |£| 2 = A, conjugated by the /i-Fourier transform. More 
specifically, we need to calculate 

(5.20) (2nh)- n h i ( f (V a V b ) dg) ((|£| 2 - (A + ^0))- 3 - (|£| 2 - (A - z0))" 3 , 1), 

Jx 

where (.,1) is the distributional pairing of a compactly supported distribution on 
with 1. The latter is, in modified polar coordinates £ = p 1 / 2 ^, \ui\ = 1, 

(5.21) \ vol(§- 1 )((p - (A + *0))- 3 - (p - (A - *0)r 3 , p^ 2 ) M . 

Since (p — (A + i0))~ 3 — (p — (A — iO)) -3 = 2ni 2~ 1 <5^ 2 ' as a distribution in p, we 
deduce that 

(5.22) @ = - vol(§"" 1 )2^ 2- 1 (n ~ 2)(n ~ 4) A^- 3 . 

Hence, 
(5.23) 

1 (n-2)(n-4) 



O = h A - n c{ J_ (V a V b ) dg) J </>{X)X%- 3 dX, c = (2tt)-" ^ ^ vol^ 1 ). 

In the general, higher dimensional (n > 8), case one simply has to keep more 
terms, one obtains more commutators, and arranges the factors similarly to how it 
was done above. Then one still needs to compute the asymptotics of all the a priori 
non-trivial terms, but at this point the only resolvent involved is the free one, and 
it is easy to see explicitly that none of the terms except the first one, namely the 
one k ept in ( 5.19 ), contribute to do, so the formal calculation that would lead to 
( 5.20| ) is indeed valid. 

So far we have only dealt with the 'averaged asymptotics', i.e. those involving 
<j)(H/X), etc., rather than with the asymptotic behavior of the C°° function a as 
A — > +oo. To analyze the latter, recall that ||i£(A + iO; h)\\ B ^ H r,s H r+2,-s^ < C/h for 
all r and for all s > 1/2 on any appropriate weighted spaces, A in a compact subset 
of (0, +00). Now, 

(5.24) R(X) = {H - X)- 1 = h 2 (h 2 H - X/h 2 )- 1 = h 2 R(X/h 2 ; h), 
so taking h = A -1 / 2 , we deduce that 

(5.25) \\R{X + M)\\„ IHr , s H r+2,- s ,<CX- 1/2 , A>1. 

Dl - r! A -l/2'- n A -l/2 > 

Such estimates follow from positive commutator estimates (such as the Mourre es- 
timate) |L7j , and require a non-trapping assumption for the semiclassical principal 
symbol fllTf - but this is just h 2 A since the potential is higher order, so the assump- 
tion is automatically satisfied. Microlocal versions of these estimates remain true. 



20 



ANDRAS VASY AND XUE-PING WANG 



For example, ( 3.16 ) is replaced, for suppF C (c, c') suppF C (c',c"), c < d < c" , 
\c\, \c"\ small, by 

(5.26) \\A(F; h)R(X + i0; h)A(F\ h)^^^,^ < C/h, 



and ( 3.17 ) is replaced, for suppxHsuppx = 0, suppF C (c, c' + e) suppF C (c', c"), 

by 

(5.27) \\ X A(F- h)R(X + i0: h)A(F; h)x\\ B{H „^ , s > } < C/h. 

In the high-energy version the right hand sides are replaced by CA~ 1//2 . The loss 
hr x corresponds to the fact that semiclassically the commutator of two pseudo- 
differential operators vanishes to one order higher than the product, see the remarks 
following the proof of Proposition |A.l| in the appendix. Similarly we can apply 
R(X + iO) iteratively and gain further decay in A, e.g. ||P(A + i0) 2 || < CA _1 on 
appropriate spaces. As in the preceeding section, below we apply these results to 
the resolvent (P Q — (A+iO)) -1 of the many-body Hamiltonian P a = aHL + (l — aHn) 
on X a x X a . 

The remark after (5.£) shows that for s > n/2, and for all ret, 

(5.28) \\K^ H/ x)\\ H r.-s <CX n '\ A>1, 

A -1 / 2 

and the same estimate holds for <5di ag with H^Zi% replaced by H~*[J* . Similar 
calculations also yield estimates corresponding to the wave front set, for example 
that for F with <£ supp(l - F), A(l - F) = A(l - F; A" 1 / 2 ) as in (|J), and for 
any r' , s' e K 

(5.29) \\A(1-F)K^ H/X} \\ <C\ n /\ A>1, 

and the same estimate holds for Jdiag- Let 

(5.30) suppF C (c',oo), c' < 0, |c'|small, 

(we can take, for example, d = —1/2), F identically 1 on a slightly smaller set, in 
particular on a neighborhood of [0, +oo). Write 

(5.31) Id = A{F) + A(l — F) + (Id -MP a /\)), 



where ipo(P a /X) is the operator with amplitude p in the quantization map (5.4). 

Since 

(5.32) 

|| (P a - (X + iQ))- 1 A(F)K i , (H/x) \\ < || (P Q - (A + iO))-M(F)||||i^ (ff/A) || 

|| (P a - (A + iQ))- l A[l - F)K m/x) \\ < || (P a - (A + iO))- 1 1| \\A(1 - F)K^ H/X) \\, 

|| (P a - (A + zO))- 1 (Id -MPJXW^h/X) || 

< \\(P a -{X + iO))-\id-i> {P a /X))\\\\K m/x) \\, 

we deduce that for s > n/2 + 1, and for all r, 

(5.33) H^aII^.-. =\\(P a -(\ + i0j)- 1 K^ H/x) \\ Hr ,- s <CX n / 4 -^ 2 , A>1. 



Now let F be such that suppP C (— oo,c'). The estimates (5.32), with a factor 

'A- 1 / 2 



of A(F) in front of (P Q — (A + iO)) , yield similar results, but in H r ' s l/2 for all 
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r , s , so we deduce that 

(5.34) \\A(F)K x \\ Hr , s , <CX n / 4 -^ 2 , A > 1. 

A-l/2 

Let 6 be a 2-cluster, and let x — Xb be a conic cut-off supported away from 
U a ^bdiag a , but near diag. Let \ be a similar conic cutoff, x identically 1 on suppx- 
Then for all a', 

(5.35) ||x(^WA)-^(^/A))||^„' <C\ n '\ A>1. 

Moreover, from ( 5.33| ) and since x(a(h)L + (1 — a)(-Mi?,) is Schwartz, we deduce 
that for all s' , 

(5.36) \\x(a(h) L + (1 - aJCJftJfl)^!!^,./ < CA"/ 4 " 1 / 2 , A > 1. 



Let v denote the right hand side of (4.7). Thus, 

(5.37) WxvWh*..' <C\ n '\ \\A(F)(l-x)v\\ H r>.*' <C\ n '\ A > 1. 
But 

#A - (^)a = (i^ - (A + iO))" 1 *; 

3 = (Pa - (A + iO))- 1 ^ + (Pa - (A + zO))- 1 ^ (Pa/A)(l - x)v 

+ (P a -(X + iO))- 1 A(F)(l-x)v 

+ (P a - (A + 4 0))-M(l - F)(l - x)v. 

Let F satisfy suppF c (— oo, e), e > sufficiently small, and apply x-A(Fq) to the 
previous equation. We can then estimate each resulting term on the right hand side 
in for all s' by CA"/ 4 ^ 1 / 2 . Namely the first three terms can be estimated by 

using the operator estimates of A(F )(P a — (\+iO))~ 1 , since it is applied to functions 
bounded in H^ll /2 by CA™/ 4 for all s'. On the other hand, the last term can be 
estimated by using the operator estimate of xA(Fo)(P a — (A-H0)) _1 ^4(l — F)(l — x). 
We therefore deduce that for all r' , s' , 

(5.39) \\ X A(F )(K x -(K b ) x ))\\ Hr ,, sl < CA"/ 4 " 1 / 2 , A > 1. 

A-l/2 

Combined with a similar result for (K a )\ — (K )\, and summing over a partition 
of unity Xb of a neighborhood of diag, we deduce a bound for x' A(Fq)K' x if x' is 
supported near the diagonal. Since the corr esponding bound is automatic away 
from the diagonal even for A(Fq)K\ by ( 5.27 ), we conclude that 

(5.40) \\A{F Q )K' X \\ H ^CA"/ 4 - 1 / 2 , A > 1. 

A-l/2 

Our proof of the smoothness of <r = (-K" A ,#diag) then first yields that |<r(A)| < 

CX (n-l)/2 

since the appropriate microlocal norms of K' x are bounded by C'A"' 4 ' 2 , 
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while those of <5 d ia g are bounded by C"A"/ 4 . Indeed, in (O), for F e Cf (R) sup- 



ported in (—00, e), identically 1 near 0, and for s > n/2 and A > 1, 
(5.41) 

|<(Id-Vo(Pa/A))i^,<W| < ||(Id-Vo(i'a/A))^|| H ...||<Jdia g ||if-..-. ^CA"/ 2 " 1 / 2 , 

\(A(F )K' x ,S dias )\ < \\A(F )k' x \\ H .,45 diag \\ H -„s < CA"/ 2 - 1 / 2 , 
|(^,(V»o(P«/A)-A(Fo)*)5 diag )| 

< \\K' x \\ H -,,-s\\(MPa/X)-A(F )*)S diae \\ H s,, KCX 1 ' 2 - 1 ' 2 . 

The identity XdR(X)/d\ = HR 2 (X) — R(X) also allows us to deduce that 
(5.42) \(Xd x ) k a(X)\ < C fc A (n " 1)/2 , 

hence a is a symbol (outside a compact set). Applying a perturbation series argu- 
ment, as for the 'averaged asymptotics' above with A replaced by A -HO, then shows 
that for all k > 0, a has an asymptotic expansion, up to A*™ - ^/ 2 , modulo symbols 
of order (n — k)/2. (Note that a priori there is only a gain of A" 1 / 2 , rather than 
A -1 , between consecutive terms of the perturbation series, due to (5.25); this does 
not appear in the trace due to the special behavior of the kernel near the diagonal: 
sharper (subconic) localization near the diagonal would give, even a priori, a better 
result.) Hence a is indeed a classical symbol, i.e. has an asymptotic expansion, 
with the top coefficient calculated above. In view of (5.23), we have thus proved 
the theorem from the introduction, which we now restate. 

Theorem 5.2. Suppose that the pair potentials V a are Schwartz (on X a ). Then 
the spectral shift function a is C°° on K \ A. Moreover, a is a symbol outside a 
compact set, and it has a full asymptotic expansion as X — > +00 : 

00 „ 

ct (a) ~ Af ~ 3 ~^> c ° = c ° E E / V - Vb d -9> 

where Cq — j^(n — 2)(n — 4)(27r)~ n vol(S n_1 ) depends only on n = dim Xq. 



6. Many-body spectral shift functions 

We now define a modified spectral shift function in the general ./V-body setting. 
Essentially the same proofs as above show its smoothness away from the thresholds 
and yield its high energy asymptotics, though now the combinatorial part becomes 
a little more complicated. We define these recursively for subsystems, starting with 
the free (i.e. N-) cluster. So, for <j> € C™(R), let 

T(X ,X ,X°)(<i>)= ( j>(H ), 

T(X o ,X a ,X a )(0) = cf ) (H a )- ]T T(X Q ,X c ,X c )(0). 

c<a, c^a 

The second equation can be rewritten as 

(6.2) 0(H a )=Y / T(X o ,X c ,X c )(0), 

c<a 

and it defines T(X , X a , X a ){<j)) recursively. Note that T(X , X a , X a )(<p) is an 
operator on X ; more precisely, it is in ^^^° ,0 (X ;C), though some of the elements 
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of C can be dropped. We write 

(6.3) T(</>)=T(X ,X 1 ,X l )(<f>). 

For example, if H is a two-body Hamiltonian, we get T{(j>) = 4>(H) — cf>(Ho), and if 
H is a three-body Hamiltonian we obtain the Buslaev-Merkurev expression, T(</>) = 
4>{H) — 4>(H ) — ^2^ a=2 (^(H a ) — 0(-^o))- We continue with a lemma. 

Lemma 6.1. Suppose that V a € x k C°°(X a ). Then for all a, T(X , X a , X a )((j>) 
is in ^5^°' (Xo,C) away from C a , i-e. on Co \ C a - In particular, if all V a are 
Schwartz, then for all a, T(Xo, X a , X a )(4>) is in i S~^^° ,oc ' (Xo,C) away from C a . 

Proof. This statement is empty for T(X , X°, X°)(<j>), since C \ C Q = 0. We 
proceed by induction, assuming that we have shown that for all c < a, c =/= a, 
T(X , X c , X c ){4>) is in ^^ ,00 {Xq, C) away from C c . Suppose that p £ Ch )reg , 
p £ C a ; in particular, b ^ a (for then C a D C b would hold). Then T(X , X c , X c )((j>) 
is in ^^^ o cc (Xo,C) near p unless p S C c , i.e. unless C b C C c , i.e. b > c. Now, 

T{X (h X\X a ){4>) = ct>{H a )- Y, T(X ,X c ,X c )(cj>) 

c<a, c^a c<b 

(6.4) " ' JL " 

- Y TiX^X^X^), 

c<a, c^a c^b 

and we have just seen that each term in the last sum is in ^ s ^°'°° (Xq,C) near p. On 
the other hand, let d be the maximal element with the property d < a and d < b. 
Note that this maximal element is unique, namely it is given by X d = span{X c : 
c < a, c < b} C X a n X b , i.e. X d = f|{X c : c < a, c < 6} D X a + X b , which is 
a collision plane since is closed under intersections. In particular, d ^ a since 
a £ b, hence the first sum in (|J) is J2 c <d T ( X o, <%" c )(0) = <X#d), so 

(6.5) T(X ,X a ,X a )(<f>) = tf>(H a )-<f>(H d )- Yl T{X a ,X c ,X c ){4>). 

c<a, c^a, c^b 

Since H a — H d = J2 c <a c<b^' ^ i s m x k C°°(X ) (resp. Schwartz) near p if the 
potentials V c are in x k C°°(X c ) (resp. Schwartz on X c ), hence the local nature of the 
construction of 4>(H a ) and <fi(H d ) (functional calculus and resolvent construction) 
yields that 4>{H a ) — 4>{H d ) is in H>s£°'°°(Xo, C) nearp, hence providing the inductive 
step. □ 

This shows, in particular, that T((j>) g \&g c °°' (Xo, C), and is hence trace class if 
k > n. Moreover, the map C£°(R) 3 <f> i— * tr(T(0)) 6 C is linear and continuous. 
We thus make the following definition. 



Definition 6.2. Let be a many-body Hamiltonian, and define T by (6.1)-( p.3[ ). 

Suppose that the potentials V a are symbols of order k > n on X a : V a e S (X a ). 
The modified spectral shift function a is defined, as a distribution on R, by 

(6.6) a(4>) = tr(T(0)), e C£°(R). 

For Schwartz potentials 14 (on X c ), the arguments presented in the previous 
sections apply, with the result that a is C°° on R \ A. Indeed, note first that the 
wave front set of kernel of T(Xq, X a , X a )((f>) is in the compressed conormal bundle 
of the diagonal (since the operator is a linear combination of the 4>{H C )), hence 
r = on it. Let ip e C~(R), and let $(t) = t/>(*)(* - (A + zQ)) -1 . We will write 
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T(Xo, X a , X a )(ip) even though ip is not smooth. We claim that near any p ^ C a , 
T(Xq, X a , X a )(i(>) is a sum of terms Tj, each of which satisfies that for some ctj < a, 
(H aj — X)Tj,Tj(H aj — A) are in (Xq,C) near p. The proof again proceeds 

by induction, the statement being empty for a = 0. So suppose that for all c < a, 
c =/= a, we have shown the claim, and suppose that p £ Cf, fleg . Now ([T^) becomes 

(6.7) T(X ,X a ,X a )(i>)=i>(H a )-i>(H d )- J2 T(X ,X C ,X C )(^). 

c<a, c^a, c^-b 

Each term in the sum on the right hand side can be written as a sum of operators 
T CJ with the desired properties by the induction hypothesis. On the other hand, 

(6.8) (H a -\)$(H a )-i>(H d ))=i>(H a )-Tl>(H d ) + { £ V c )iP(H d ). 

c<a, c^-d 

Since c < b, together with c < a, would imply c < d, we deduce that each V c is 
Schwartz near p, hence the last term is in *f/g£°' 00 (Xo, C) near p, while the same 
statement for ip(H a ) — ip(H d ) has already been demonstrated above. Proceeding in 
the three-body setting, i.e. using the propagation of singularities, shows that the 
wave front set of the kernel of T(ip) is in r < tq < 0, and then we deduce that a is 
continuous. Iterating the argument gives that a is C°° . 

The pseudodifferential functional calculus shows, as in the previous section, that 
for0eC c °°(A), 

OO 

(6.9) ^./A))-^"-^-, 

3=0 

with a,j = for all odd j. The symbol estimates for the C°° function a itself also pro- 
ceed as before, so \(Ad\) k a(A)\ < CfcA^"" 1 ^ 2 . A perturbation series argument again 
yields a full asymptotic expansion. Since R a (X) — Rq(X) — R a (A)(J2 c <a ^c)-Ro(A), 
A ^ K, we deduce that the leading terms ao and a\ vanish. If all potentials are 
pair potentials, i.e. if V c = for all clusters c that are not N — 1-clusters, one 
can show that a,j — for j < 2N — 1 and ci2N = Cti.nC-(V) where C(V) = 
J2 a J Vai ' ' ' Vcr N -idx, the sum is taken over all sequences a — (a%, • • • , <tn-x) of 
(N — 1) -cluster decompositions Oj such that U^T^aj = a max and c n .N is a contant 
depending only on n and N which can be calculated by the method of Section [|. 
We have thus proved the following result. 

Theorem 6.3. Suppose that for all a, V a is a Schwartz function on X a . The 
modified spectral shift function, a, defined as a distribution by ( |6.6| ), is C°° on 
K\A. In addition, a is a symbol outside a compact set, and it has a full asymptotic 
expansion: 

oo 

<j(\) ~ ^2 X^~ 2 ~ j Cj, A^+oo. 
j=o 

Here we do not assume that all potentials are pair potentials. 

Appendix A. Sketch of relevant positive commutator estimates 

To illustrate the proof of the propagation of singularities by positive commutator 
estimates in pi) , we sketch the proof of a simpler version here which still suffices 
for the purposes of the present paper. We state it as a wave front set estimate, 
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though as we see below, as usual with positive commutator estimates, it actually 
amounts to a microlocal energy estimate. 

First we introduce some notation. For an cluster a (possibly a = 0!) let 

(A.l) x a = \wa\~ 1 , y a = w a /\w a \, z a = w a / \ w a | ; 

these are local coordinates on the radial compactification Xq near X a . Let 

(A. 2) T = j |-, T a = j j , 

so r = r a at X a . We also write x = In addition, it is sometimes convenient, 

for AT b D X a , i.e. X b C X a , to decompose w a € X a as (w ab ,w b ) € X b ®(X a eX h ), 
and write (z ab ,z b ) accordingly. We also let sc Hg be the rescaled Hamilton vector 
field of Ax c , so 

(A.3) ^ H c 2 ^c. d 

y \Wc\ 

this should be regarded as a vector field on T*X C which extends to a smooth vector 
field on SC T*X C tangent to the boundary SC T<X X c , and hence can be considered as 
a vector field on SC J£ X c . Recall also that the part of the characteristic variety 
corresponding to the bound states of H b is 

(A.4) E 6 (A) = {C = (2/6,6) G sc T* c X b : A - |6| 2 G spe Cpp H b } c sc T* c X b . 
For C a C Cb, we write the projection 

(A.5) 7r ba : £ b (A) n sc T^X b -» sc T^X a c E(A); 

this is the restriction of n ba : sc T£ ai ^X b -> sc T£ a t X a to £&(A). For A e 

^g,, 00 ' (X,C), the operator wave front set WFg c (A) was defined in |3(| as a subset 
of the compressed cotangent bundle SC T*X. Namely, let p denote the projection 
[X;C] x X* — > SC T*X given by the composition of the blow-down map and the 
projection tt : SC T*X -> sc f*X. Then C G s T*l\WF' Sc (i) means that in a neigh- 
borhood of p~ 1 ({C}), the amplitude a defining A as in ( |3.6[ ) vanishes to infinite 
order. This notion is independent of the choice of quantization. 
The main technical result is thus the following. 

Proposition A.l. pl| , Weaker version of Proposition 7.1] Suppose that H is a 
many-body Hamiltonian. Let u € C~°°(X), X ^ A\. Let y a G C a . reg , f such that 
A — f 2 ^ Ai. Suppose that for all C = {y a Aa) G SC T£ X a with r(Q = f, we have 
C ^ WF sc({H — X)u). Then there exist So > and C > 0, depending only on 
wii/i f/ie following property. For all Sq G (0, <5o) sucft i/iai z/ 

(A.6) VC s.t. 1^(0 ~ Val < CSo and f + S /3 < r(0 < f + «J C £ WF Sc ( U ), 

i/ien C ^ WF Sc («). 

In fact, there exists C > so that the following holds. For any r, s € M, r' > — r 7 
s' > — s, depending only on £ 7 f/iere exists C\ > suc/i i/iai /or u as above, 

(A.7) 

ll^(XC/'(yo)X(-«o/3,5o/3)(^-T))w||j?^ < Ci(||A(xc/(2/o)X(5o/3^ )( T -' r )) w llff''. s 

+ \\Mxu{ya)X(-s a /3,s )( T ~ f ))( H - X)u\\ B r,.+i + \\u\\ B r>,.'), 
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where U is the ball |y (C) ~ Va\ < CS , U' the ball |y (C) ~Va\ < C'S , xu(y), resp. 
Xi(t — f) denote a smoothed characteristic function ofU y , resp. the interval I T -f, 



and A denotes quantization as in (3.9) 



Remark A. 2. Since £ £ WFg c ((i? — \)u), by elliptic regularity we deduce that 
C 4- WFsc(u) for £ £ S(A), i.e. for ( not in the A characteristic set (energy shell). 



The estimate of (A. 7) implies (3.17) directly. Indeed, if x is such that suppxu 1 H 
suppx = 0, consider u = R(X + iO)v, v = A(x(yo)X(-5 /i,S /4) ~ f))f, f G H r '' s '. 
Then the second and third terms are directly bounded in terms of v, while the first 
term is bounded in terms of / due to the boundedness of -<4.(X(*o/3,*o) ( r ~ ^))-^(^ + 
iO)A(x(-deita /4,5 /4) ( T — t)) between any two weighted Sobolev spaces. 

Proof. (Sketch, see Proof of Proposition 7.1] for complete details.) We give 
the full commutator construction at the symbol level, and indicate why it gives rise 
to a microlocally positive commutator. In fact, the commutator will be positive in 
part of phase space, negative (or not necessarily positive) in another part of phase 
space. The propagation of singularities estimates, which should be thought of as 
microlocal energy estimates, work by estimating u in the former region in terms of 
u in the latter region and (H — X)u in the union of both regions. 

Employing an iterative argument, we may assume that for all C — (j/aiCa) with 
r(C) = r, C i WFgf (u), and we need to show that C £ WFg^ +1/2 (u). (We can start 
the induction with an I such that u £ H*' l (X).) 

For points C = (ift,&) in sc T^X b , 

(A.8) - sc H b g r(C)=2(\^\ 2 -r(Cf)>0. 
For C € E 6 (A) c sc T* b X b , A - |^| 2 =^6 spec pp (H b ), so 
(A.9) - sc H b g T(C) = 2(\-T((f-e ). 

We define 
(A.IO) 

co = \ mi{- sc H b g T(( b ) : 3£ a such that r((y B ,€„)) = f, ( b € it£{(y a ,£ a )), C b D C a }. 

Due to ( |A.9| ), and due to A — f 2 £ A l5 we deduce that cq > 0. Thus, there exists 
Si > such that for all clusters b with C b D C a , and for all £ € E&(A) that satisfies 
\y a (0 ~ Va\ < \ z a{()\ < Si, |r(C) - f | < Si, we deduce that 

(A.ll) - sc H b T(C) > 3c /2. 

Our positive commutator estimates will arise by considering functions 

(A.12) 4> = f -T+P(\z a \ 2 + \y a - y a \ 2 ), 

e 

where (3 > will be fixed later and e > is arbitrary as long as it is sufficiently small. 
Note that for all b with C b D C a , sc H b \z a \ 2 = 4z ab ■ under the decomposition 

= (£a,£&)j so sc H b \z a \ 2 < Ci\z a \ on E b (A), and similarly, possibly by increasing 
Ci, \ sc H b g (y a -y a ) 2 \<Ci\y a -y a \. 

Now suppose that 

(A.13) 4> < 2e, f-r > -2e. 
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Then we conclude that 

(A.14) |f - r| < 2e, \z a \ < 2e/y/0, \y a ~Va\< 2e/y/0. 

Let = (c /4Ci) 2 . For e > small, ( |A.13| ) th us im plies that \y a {() - y a \ < S ± , 
\ z a(C)\ < Si, \t(Q — f | < 8%, so we deduce from ( A.ll ) that 

(A.15) sc H b g cf> > - sc H b g T - > c , where (3 = (c /4Ci) 2 . 

The positive commutator estimate then arises by considering the following sym- 
bol q and quantizing it as in (B.E). Let xo € C°°(R) be equal to on (— oo,0] and 
Xo(t) = exp(-l/t) for * > 0. Thus, x' (t) = t~ 2 xo(t), t > 0, and X 'o(t) = 0, t < 0. 
Letxi eC°°(R) beOon (-oo,0], Ion [1, oo), with X 'i > OandxiO) = exp(-l/i) on 
some small interval (0, to), to > 0. Furthermore, for Ao > large, to be determined, 
let 

(A.16) q = xo(Ao 1 (2 - <P/e)) X i((f - r)/e + 2). 

Thus, q(C) = xo(2/Aq) > 0, and on suppg we have 
(A.17) < 2e and f - t > -2e, 

which is ( A.13Q , so suppg is a subset of (A.14). We also see that as e decreases, so 
does suppg = suppg e , in fact, if < e' < e then q c > on suppqv. Note that by 
reducing e, we can make q supported in an arbitrary small neighborhood of y a and 
f a - 

Let t/j € C^°(K) be identically 1 near and supported close to 0. We also define 
(A.18) q = $(x)q. 



Let A be the operator given by (3.£) with q in place of q. Note that this includes 
a spectral cutoff in the definition of A. 

The commutator i[Ax c , A] is given to top order by sc H g q. This is the com- 
mutator that gives microfocal positivity on the L 2 eigenspace of H c , see e.g. the 
Froese-Herbst proof of the Mourre estimate We proceed to esimate sc H g q di- 
rectly. 

Thus, 

(A.19) 



= -A^e-Wo(Ao 1 (2 - <f>/e)) X i((f - r)/e + 2f c H^ 
- e-\o(A^(2 - 0/e))xi((f - r)/e + 2) sc H c q 



Then 

(A.20) sc H c g q = -bl + e c 
with 

(A.21) b 2 c = A^e' 1 X ' Q (A^(2 - 0/e)) Xl ((f - r)/e + 2) SC H^. 
Hence, with 

(A.22) b 2 = coA^e- 1 X 'o(Ao 1 (2 - 0/e)) X i((r - r)/e + 2), 
we deduce that 

(A.23) sc H c g q < -b 2 + e c . 
Moreover, 

(A.24) b 2 > (c Ao/16)q 
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since <f> > r — r > — 2e on supp q, so 

x'o(V(2-^A))=^(2-0/e)- 2 X o(V(2-0A)) 

(A.25) 

>(^g/16)xo(V(2-0A)). 
On the other hand, e c is supported where 

(A.26) -2e < f - r < -e, |y a - y a \, \z a \ < 2e/y/p. 

By our assumption, this reg ion is disjoint from WFs c (u), if we choose e > suffi- 



ciently small. Moreover, by ( A.f 7 ), for e > sufficiently small, we deduce from the 
inductive hypothesis that supp q (hence supp b) is disjoint from WF*;'(«) n S(A). 
Let B G Vs™ fi {X,C) be a quantization of bq 1 / 2 as in |P|). Suppose that M > 



and e' > 0. By choosing A large, depending on M, e', (using ( A.25| )), one can 



derive a positive commutator estimate from (A.23) using the many-body pseudo- 
differential calculus, see |3l], End of proof of Proposition 7.1] for details. Apart 
from technical details it essentially corresponds to using the Mourre estimate and 
the functional calculus microlocally, namely that when localized in phase space in 
the region of interest, the commutator of a quantization of <fi is positive. We deduce 
that there exists 5' > 0, such that for ip e C^°(M) is supported in (A — 5', A + §'), 
= 1 near A, E G *s c ° o,0 (X,C), F G 9^° A {X,C) with WF' Sc (E),WF' Sc (F) in a 
small neighborhood of S(A), 

(A.27) WF' Sc (£) c suppe, WF^ C (F) c suppg, 

such that 
(A.28) 

i^{H) x - 1/2 [A*A, H] x - 1/2 iP{H) - Mip(H)A*Aip(H) > (2 - 2e')ij)(H)B* Bip(H) + E + F. 
By supp e we mean the support of the function xoi^Q 1 ^ ~ < / ) / e ))xi((^ — T )/ e + 2), 



which is independent of c in (A. 19). Here F is the error term, it has first order 
decay, hence it is 'negligible'. On the other hand, E has the same order as B*B, 
and it is negative (i.e. has the opposite sign of B* B) in part of the phase space. 
As mentioned above, positive commutator estimates for approximate solutions u, 
i.e. (H — X)u microlocally Schwartz, work by estimating ||S^(i?)u|| 2 in terms of 
{u, Eu) (plus error terms), i.e. u is estimated on supp b by its estimate on suppe. 

One can now use M, chosen sufficiently large, to deal with arbitrary weights 
x -l-i/2^ standard commutator and regularization argument then proves that 
x~ l ~ 1 / 2 Bu G L 2 {Xq), which in turn finishes the proof. We refer to |3l], Proposi- 
tion 7.1] for details. 

Instead of following this route, we prove the corresponding resolvent estimate. 
So suppose that 

(A.29) uf = (if — (A + it))- 1 !, t > 0, 

and WFg c (/) is disjoint from the region of interest, and it is, say, in r < t , t > 
sufficiently small, so that uf converges to (if — (A + iO)) -1 / as t — > in sufficiently 
large weighted Sobolev spaces. As above, assume that uf is uniformly bounded 
in the region of interest in H* ,l (Xo); we want to prove that it is also uniformly 
bounded in if *'' +1 / 2 (A ). For ip G C™(R; [0,1]) supported sufficiently close to A, 
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with Ai = Aif)(H)x 
(A.30) 



-1-1 



Bi 



-1-1/2 



Btp(H), we deduce from (|A.28|) that 



/./- ■'~[A^A l ,H}x l+1/2 > x l+1 l 2 {{2-2e')B* l B l +E l +F l )x l+1 l 2 , e' > 0, 



(X,C), Fi £ * s °°' 2l (X,C), with similar properties as in flA.27p. 



1+1/2 
-oo,-2/-l 

Since 

(A.31) (u+,i[A*Ai,H]u+) = -2Im(u+ A* A t {H - (\ + it))u+) -2t\\Aiu+f. 

we conclude that 
(A.32) 

\\Biu 



+ l|2 



2t\\A lU t\\ 2 < \{u+,E lU t)\ + \(u+,F l u+)\+2\(u+,A^A l (H-(X + it))u+)\. 



Since t > 0, the second term on the left hand side can be dropped. Since — > u+ 
in H"- l '(X Q ) for I' < -1/2, we conclude that for / S (-1,-1/2) the right hand 
side stays bounded as t — > 0, for uf is uniformly bounded in H°' l+1 / 2 (Xo) on 
WFg c (£; ; ) and it is uniformly bounded in H a > l {X a ) on WFg c (F/). Thus, B t u+ is 
uniformly bounded in L 2 (Xo), and as — > u+ in (Xq), we conclude that 
G L 2 (X ). 

□ 



The semiclassical version of the estimate (A. 7) is 
(A.33) 

\\A(xu>(yo)X(-s /3M/3)(T - f))u\\ H r, s < Ci(\\A(xu(yo)X(s /3,s )( T ~ f )) u \\H r h s 

+ h- 1 \\A{xu(yo)X(~8 /3.8 )( T - f ))( H - A)u||„,, s + i + Hl^'.a')- 



The proof of this proceeds just as above. Equation ( A. 28 ) is replaced by 

i4>(H)x- 1/2 [A* A,H{h)]x~ 1/2 %ij{H) - hMip(H(h))A*Aip(H(h)) 
> (2 - 2e')hijj(H(h))B*Bip(H(h)) + hE + h 2 F, 



(A.34) 



i.e. the principal terms B* B and E have an extra factor of h (since the semiclassical 
calculus is commutative to top order in h), and the error term F has a gain of h. 
Note that H(h) = A + h 2 V shows that semiclassically V is two orders lower in h 
that A, which in fact significantly simplifies the argument that turns (A. 23) into 
(A. 28) (for sufficiently small h). Then (A.32) becomes, after dropping the second 
term on the left hand side and multiplying through by h^ 1 , 

(A.35) 

||£^+|| 2 < \(u+,E lU +)\ + h\(ut,F lU +)\ + 2h- 1 \(u+ 1 A* l Ai(H - (A + it))ut)\, 



and then one can finish the proof as before, using that Uj is bounded by Ch 
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